Astronomy <fe Astrophysics manuscript no. nicmos2'08'nov'2004 


February 2, 2008 


(DOI: will be inserted by hand later) 





The new infra red view of evolved stars in I Zw 18 * 

G. Ostlin 1 and M. Mouhcine 2 ' 3 



o 
o 

(N 

> 
O 



> 

(N 



O 

6 



X 



Stockholm Observatory, AlbaNova University Center, SE-106 91 Stockholm, Sweden 
e-mail: ostlin@astro.su.se 

2 School of Physics and Astronomy, University of Nottingham, Nottingham NG7 2RD 

3 Observatoire Astronomique de Strasbourg (UMR 7550), 11, rue de l'Universite, 67000 Strasbourg, France 
Received 28/04/2004; accepted 31/10/2004 

Abstract. We report results from near-infrared imaging of I Zw 18, the most metal-poor galaxy in the local 
universe with NICMOS on board the Hubble Space Telescope (HST). Observations were obtained in the broad 
F205W filter (~ K) and in the medium-wide bands F171M and F180M, the latter which cover strong molecular 
absorption features in cool stars. The new data, together with previously obtained HST/NICMOS images in the 
F110W (~ J) and F160W (~ H) filters, provide a census of the cool stellar population in this chemically unevolved 
galaxy. We find that stars as old as ~ 1 Gyr are required to explain the observed colour-magnitude diagrams. 
Combining broad and medium-wide band photometry, we have classified the observed stars. The observed stellar 
populations, to the depth of our data, are dominated by luminous red supergiant stars. However, an intermediate- 
age component is also present. We have identified carbon star candidates, and show that they dominate the 
stellar content for intermediate mass and age. We show that low metallicity intermediate-age stars have, for a 
given colour, a F110W (~ J-band) excess compared to local carbon stars. The carbon stars identified in IZwl8 
are the most distant (more than 10 Mpc) resolved examples yet discovered. 
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1. Introduction 

There has been a strong interest in blue compact galax- 
, ies ever since the pioneering study of IZw 18 by Searle & 
Sargent (1972J). The global properties of IZw 18, i.e., ex- 
| tremely low nebular oxygen abundance, 12 + log(0/H) = 
7.18 (Izotov and Thuan H999I just a few percent of the 
solar abundance), very blue colours, and a large Hi reser- 
voir, have been interpreted as evidence that this galaxy 
may be genuinely young, presently forming its very first 
generation of stars. If its star formation history really is 
restricted to the last few 10 Myr, the existence of I Zw 18 
demonstrates that the formation of galaxies may be con- 
tinuing also at the present cosmic epoch. Speaking against 
the youth hypothesis is the discovery of a population of 
faint red stars in I Zw 18, both in the optical (Aloisi et al. 
irgfjgjl and in the near-infrared (Ostlin 2000). These stars 
have probable ages in the range from 100 Myr to a few 
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Gyr, and most of these are expected to be stars on the 
asymptotic giant branch (hereafter AGB). 

The upper AGB is populated with oxygen-rich and 
carbon-rich stars (i.e., AGB stars with C/O ratio smaller 
or larger than unity by number, respectively). It has been 
established observationally and theoretically that carbon 
star formation is more efficient in low metallicity envi- 
ronments. The number ratio of carbon stars to M stars 
anti-correlates with the metallicity (Cook et al. 1986, 
Groenewegen Unnni Mouhcine & Lancon 2003). IZw 18 
has the lowest ever measured ionised gas phase metallic- 
ity in a galaxy (Kunth & Ostlin 2QQ0J. Thus, if there is 
an intermediate-age stellar population present in I Zw 18, 
one might expect it to be dominated by carbon stars. In 
the same context, the carbon-star content of I Zw 18 could 
provide insights into stellar evolution at low metallicity. 

AGB star spectra are dominated by broad molecular 
absorption bands. Spectra of oxygen-rich stars are domi- 
nated by metal oxide bands such as TiO, VO, and H 2 0. 
Optical and near-IR spectra of carbon stars have little 
in common with those of oxygen-rich AGB stars apart 
from the CO absorption in K-band. They are dominated 
by bands of CN and the Phillips & Ballik- Ramsay C2 sys- 
tems (Barnbaum et al. 1996, Joyce 1998). Because of their 
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high luminosity and distinctive spectra, carbon stars can 
in principle be identified out to great distances. It has 
been argued that carbon stars have, in general, redder 
broad band colours than oxygen-rich stars. This statement 
is not entirely true. For AGB stars with low mass loss rate, 
i.e., stars with optically thin circumstellar shells, carbon 
stars display a wide range of colours, similar to the colour 
range spanned by M-stars (Alvarez et al. 2000). Due to 
CN and C2 absorption in the J-band, the J — K colour 
distribution of carbon stars extends to redder colours than 
M-stars, i.e., J - K ~ 2.5 (Loup et al. 1998). Applying 
a sharp cutoff to separate carbon stars from oxygen-rich 
stars may miss warm carbon stars, and misclassify ex- 
treme oxygen-rich optically selected or dust-enshrouded 
stars as carbon stars. Thus, optical/near- IR colours are 
insufficient, by themselves, to discriminate carbon stars 
from oxygen-rich stars. 

A more successful strategy to find carbon stars is to 
make use of spectral features that uniquely identify such 
stars. Different indices, corresponding to different spec- 
tral features, have been proposed in the literature (Wing 
1967, Aaronson et al. 1982). The NICMOS instrument on 
the Hubble Space Telescope (HST) has a set of medium- 
wide filters. One of these, F180M in the NIC2 camera, 
coincides with a C2 band head at 1.77/im (Ballik- Ramsay 
series, Aw = 0). Together with the adjacent continuum 
filter F171M, it can be used to probe the strength of the 
C 2 band head (Alvarez et al. 2000). 

In this paper, we report the results from 
HST/NICMOS medium-wide (F171M and F180M) 
and broad (F205W, ~ K) band imaging of I Zw 18, which 
we use to investigate the nature of red stars and the 
upper AGB population. The observations and the data 
reductions are described in Sect .[21 In addition we use 
previously obtained NICMOS data in the F110W (~ 
J) and F160W (~ H) filters. In Sect.[3 we discuss the 
surface photometry. The nature of the red and cool stellar 
content of I Zw 18 identified using colour-magnitude and 
colour-colour diagrams is discussed in Sect.0J Conclusions 
are summarised in Sect.|S] 

Table 1. Log of observations for GO programme 7880, 
with zero points given in the Vegamag system. 



Filter 


Total exposure time 


Zero Point 




object 


sky 






(8) 


00 




F171M 


4629.47 


3086.32 


19.956 


F180M 


6172.63 


4115.09 


19.844 


F205W 


4629.47 


3086.32 


21.960 



2. Observations and point source photometry 

The observations were acquired using the NIC2 camera of 
HST/NICMOS under HST general observer programmes 



Fig. 1. RGB composite of F205W (red), F160W (green) 
and F110W (blue) showing a field of 14" x 14". The orien- 
tation is indicated by the green compass in the lower right 
corner (north is up right, east is up left). The length of the 
bars is 2". Two of the diffuse sources with strong excess 
in the F205W filter are indicated by thin white arrows. 



7461 and 7880. For a description of the NICMOS instru- 
ment, see Thompson et al. I)1997|l . The observations for 
programme 7461 consisted of a total of 2560 s exposures in 
the F110W filter and 5120 s in F160W, and are described 
in Ostlin ( 2000 ) together with details on the reduction and 
data analysis. 

The observations for programme 7880 were carried out 
during 10 orbits and are summarised in Tabled Due to the 
long wavelengths of the filters, the observations were ob- 
tained by using "chopping" , i.e. separate sky frames were 
obtained in between each exposure on target. The obser- 
vations employed the SPIRAL-DITH-CHOP pattern in the 
MULTIACCUM mode and the exposure sampling sequence 
used was MIF512. Hence, for each orbit, five 512-second 
exposures were obtained, three on the primary target and 
two on adjacent sky positions. In order to improve the 
spatial sampling and photometric accuracy, small (1.5") 
offsets were made in between the individual exposures for 
both the target and sky positions. Some images suffered 
from the so called pedestal effect, which was removed us- 
ing the Pedestal Estimation and Quadrant Equalization 
Software, developed by R.P. van der Marel. 

Before combining the target images, the background 
was subtracted using the sky frames. Best results were 
obtained by first combining all sky frames to a master 
sky which was then subtracted. The sky-subtracted target 
frames were finally combined and rebinned to a finer pixel 
grid using the drizzle package (Fruchter & Hook, 1998). 



Ostlin & Mouhcine: Carbon stars in I Zw 18 



3 



Point source photometry was performed using the 
DAOPHOT package (Stetson et al. 11990(1 . Aperture cor- 
rections to convert our daophot magnitudes to total 
NICMOS magnitudes 1 were determined by performing 
photometry of synthetic NICMOS PSFs generated with 
the Tiny Tim software (Krist & Hook 11557} ! . 

The observations were acquired at a different space- 
craft roll angle compared to programme 7461, and the 
astrometry of point sources for the two different programs 
were combined with the xy2rd task in STSDAS. In match- 
ing the photometry catalogues in different filters, a tol- 
erance of 0".0375 was applied for the centre coordinates. 
Compared to Ostlin (2000), we detect somewhat fewer 
stars due to the narrower band passes of the F171M and 
F180M filters, and higher background in F205W. The zero 
points are given in Table ^ and are in the Vegamag sys- 
tem. We will denote magnitudes in the different filters 
simply by their name, e.g. F205W. 

Table 2. Model and observed colours of extended emission 
features. The model colours have been calculated with the 
Zackrisson et al. (2001) code and assume a burst with age 
3 Myr, a standard Salpeter IMF with a mass range of 0.08 
to 120 M , and gaseous and stellar metallicities of Z stai s = 
0.001 and Z gas = 0.0004. Note that the effect of nebular 
emission on broad band colours sensitively depends on the 
assumed metallicity and the filter set used. The model 
colours presented here are only valid for I Zw 18 and the 
NIC2 filter set. 



Colour 




Model 
stars 


stars+gas 


Filament 


Blob 


F110W- 


-F160W 


-0.05 


0.23 


0.5 


0.3 


F110W- 


-F205W 


-0.06 


1.13 


1.3 


0.9 


F160W- 


-F205W 


-0.01 


0.90 


0.7 


0.6 



3. Extended and nebular emission 

The galaxy consists of two main star-forming complexes: 
The bright north-west (NW) and the more diffuse south- 
east (SE) regions, which are embedded in a faint irregular 
envelope. When analysing the F205W (~K) vs the F110W 
(~J) and F160W (~H) images some features become ev- 
ident (see Fig.Q. To the west of the NW region there is 
a filament, and approximately 1" to the east of the centre 
of the NW complex there is a "blob" with size ~ 0.5" . 
Both features are bright in the F205W filter but rather 
faint in both F110W and F160W. Their colours are given 
in Table 

These features, and other larger filaments have also 
been seen in optical Ha images (Dufour and Hester HMH 
Hunter and Thronson [1995| Ostlin et al. 119961 Cannon et 

1 NICMOS total magnitudes are defined as the magni- 
tude within an 0.5" radius aperture plus —0.156 mag.; see 
www.stsci.edu /instruments /nicmos 



al. I2002|l . The filament corresponds to the regions NW D4 
and NW D5 in the nomenclature of Cannon et al. ( 2002 ) , 
and the blob to NW 1 , all of which have very small red- 
dening. 

The effect of nebular emission on the NICMOS colours 
has been computed for a stellar population including neb- 
ular emission using the models of Zackrisson et al. (2001). 
The realistic treatment of the nebular emission and the use 
of the HST+NIC2 throughput curves make these models 
suitable for our current purposes (see Ostlin et al. 2003). 
We assumed a Salpeter initial mass function with mass 
range 0.08 to 120 M Q , Z stars = 0.001, Z gas = 0.0004, a 
hydrogen density of 100 cm~ 3 , a filling factor of 0.1 and a 
gas covering factor of unity. Since the near-IR spectrum is 
completely dominated by nebular emission (line and con- 
tinuum) for ages less than a few Myr (see e.g. Bergvall 
& Ostlin 2002) the composite stellar+gaseous spectrum is 
a good approximation to a pure nebular spectrum. The 
model colours with and without gas emission for an age 
of < 3 Myr and constant star formation rate are given in 
Table Given the photometric uncertainties, the colours 
of the filament and the blob are consistent with pure nebu- 
lar emission. The high Ha surface brightness suggest that 
the blob is ionised from inside by a cluster of young stars, 
or alternatively could be a supernova remnant. The ef- 
fects of nebular emission on the colours for IZwl8 have 
also been discussed in Izotov et al. (2001), Papaderos et. 
al. (2002) and Hunt et al. (2003). 

The possible contamination of point source photome- 
try by patchy nebular emission has been assessed by com- 
paring the spatial locations of our stars with the Ha maps 
in Cannon et al. (2002). None of our carbon star can- 
didates (see Sect. 4) with available F171M and F180M 
photometry are located near Ha sources or regions with 
high Ha equivalent width. Of the AGB star candidates 
without F171M and F180M photometry that are clas- 
sified solely based on the F110W, F160W and F205W 
colours, about five are located near, but not coincident 
with, Ha sources in the SE complex. Hence, we do not 
expect our point source photometry to be contaminated 
by nebular emission. As discussed in Sect. 4.1, if con- 
tamination by nebular emission was an issue, the effect 
would be that the intermediate-age population would be 
older than it now appears. Moreover, nebular gas with 
the metallicity of I Zw 18 has medium wide band colours 
F171M-F160W > -0.2, and F180M-F171M < 0, in- 
compatible with the colours of our carbon star candidates 
that we discuss below. 

As described in Sect. 0] carbon star candidates in 
IZw 18 are characterised by their location in the F171M- 
F160W vs F180M-F171M colour-colour diagram. To test 
whether we could see any evidence for carbon stars with 
luminosity below our limiting magnitude for PSF photom- 
etry, we investigated the distribution of individual pixels 
in these colours. This was done after binning the image 
to the original pixel scale of NIC2 (0.075") and applying 
a lower flux limit of 26.6 magnitudes per binned pixel in 
F171M and F180M, in order to omit too uncertain pix- 
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Fig. 2. Observed colour-magnitude diagrams: F205W vs. F110W-F205W (left) and F205W vs. F160W-F205W (right). 
Isoclirones for a metallicity of Z — 0.0004, and a distance modulus of [m — M) = 30.5 are overplotted. Each isochrone 
is labelled with its age. 



els. Pixels containing stars already described by the single 
star PSF photometry were also omitted. In the crowded 
NW region, most pixels are contaminated with light from 
bright stars, but in the SE region we see tentative evi- 
dence for carbon stars below our point source detection 
limit (see Fig. Hand Sect. l4.5|l . 

4. Discussion: the nature of infrared-detected 
stars in I Zw 18 

In the following subsections we will discuss the results of 
the resolved star photometry. By comparing with model 
isochrones and empirical data of evolved late type stars, 
we constrain the nature of of the faint red stars in I Zw 18. 

4.1. Colour-magnitude diagrams and model isochrones 

In Figure El we show the F205W vs. F110W-F205W, 
and F205W vs. F160W-F205W colour-magnitude di- 
agrams (hereafter CMDs). Similar to what was found 
by Ostlin (2000), the CMDs display two main features: 
The first feature consists of bright (F205W < 22) stars 
with relatively blue colours (F110W - F205W < 1.2, 
F160W - F205W <, 0.3). The second feature con- 
sists of faint (F205W 22) stars with red colours 
(F110W - F205W > 1 and F160W - F205W > 0.2). 
This suggests that the star formation history of I Zw 18 
has involved two separate (in time) events or has been 
continuous over a longer time. 

A means for interpreting the observed CMDs is pro- 
vided by isochrones calculated from the stellar popula- 



tion models by Mouhcine (2002) and Mouhcine & Lancon 
(2002), which include the evolution of stars through the 
upper AGB phase. The theoretical predictions of these 
stellar population models have been extensively and suc- 
cessfully compared to a variety of observational con- 
straints. In order to constrain the stellar content of I Zw 18, 
a set of isochrones with a metallicity of Z — 0.02 Zq were 
calculated. For a detailed discussion of the models, their 
ingredients and limitations, the reader is referred to the 
papers cited above. An advantage of these models is the 
inclusion of carbon stars in a self-consistent manner, ac- 
counting for the dependence of mass loss and the sensi- 
tivity of carbon star formation to the stellar mass and 
abundance. This is of crucial relevance here as our aim 
is to identify carbon star candidates. To transform the 
theoretical luminosity-effective temperature diagram to a 
colour-magnitude diagram, we used the library of Lancon 
& Mouhcine (2002) for AGB stars, and the compilation 
of Lejeune et al. (1998) for all other types of stars. The 
former library contains spectra of both oxygen-rich stars 
and carbon stars. The isochrones have been calculated for 
the HST/NICMOS photometric system, with the aid of 
the NIC2 throughput curves. 

To compare the predicted isochrones to the observed 
colour-magnitude diagram, a distance to I Zw 18 is needed. 
This is still a matter of debate, varying according to dif- 
ferent authors from 10 Mpc to 20 Mpc. Hereafter, we will 
adopt a distance modulus of (rn — M) = 30.5 for I Zw 18 
(see Ostlin 2000 for a discussion of the choice of this dis- 
tance). Internal reddening by dust has been found to be 
very low (Cannon et al. l2002|l . and has been neglected. As 
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a word of caution we note that even if the average dust 
content is low, there are regions with a locally higher ex- 
tinction (Ay < 0.5 mag). These are however not spatially 
coincident with the sources of primary interest in this pa- 
per and moreover, the effect on the infrared colour indices 
under study here are marginal (for instance Ay = 0.4 mag 
corresponds to E(F160W - F205W) ~ 0.02). 

Figure^] shows a set of isochrones, with ages noted 
in each panel, superposed on the observed CMDs. The 
isochrones show that once AGB stars appear in a stellar 
population, at ages older than ~ 100 Myr, a red tail de- 
velops in the CMD and extends to quite extreme colours. 
As a consequence of high carbon star formation efficiency 
at low metallicities, AGB stars turn into carbon stars al- 
most immediately after they enter the red tail. The extent 
of the red tail, and particularly the carbon star sequence, 
depends on the calibration(s) used to transform the fun- 
damental stellar properties of the AGB sub-population 
stellar to observed quantities. Using the stellar library 
of Lancon & Mouhcine (2002), the portion of isochrones 
with colour redder than F110W-F205W ~ 1.4 and 
F160W — F205W ~ 0.3 are occupied by carbon stars. 

For a distance modulus of (m — M) = 30.5, the ob- 
served bright plume with its relatively blue colours, is con- 
sistent with a young stellar population in the age range 
10-30 Myr, corresponding to the age of turn-off stars with 
masses in the range of 21-10 M©. The fainter end of the 
this plume corresponds to a slightly older population, with 
age ~ 50 Myr and a turn-off mass ~ 6M Q . To account for 
the faint and red sequence, stars in the age interval 100- 
800 Myr are required. Changing the distance to the galaxy 
within a reasonable range can not alleviate the need for 
stars older than 100 Myr. All the isochrones of stellar 
populations younger than a 100 Myr predict a plume of 
stars with colours bluer than (F110W - F205W) « 1 and 
(F160W - F205W) « 0.2. To occupy regions in the CMDs 
with redder colours than these limits requires the pres- 
ence of intermediate-age stars with turn off masses in the 
range 5-2 M@ evolving through the thermally pulsating 
AGB phase. This suggests that IZwl8 started forming 
stars at least 0.5 Gyr prior to the current star formation 
event. Stars that are located in the CMD regions occupied 
by faint and red stars are, according to the isochrones, 
mostly carbon stars. 

A comparison with the theoretical isochrones shows 
that while the intermediate-age stars in I Zw 18 span a sim- 
ilar range in F160W-F205W to the models their F110W- 
F205W colours are on average bluer. This cannot be 
fully explained by the effect of a low initial metallicity of 
intermediate-age progenitors on the effective temperature, 
as this would make both colours bluer. 

One possibility is that the colours of intermediate-age 
stars have been contaminated by nebular emission in the 
neighbourhood. The effect of adding nebular emission is 
to make the F110W-F160W colour appear bluer than 
its intrinsic value, whereas the effect is much smaller for 
F160W-F205W since for this colour index, nebular gas 
and AGB stars have similar colours. However, as discussed 
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Fig. 3. F110W-F160W vs. F160W-F205W diagram for 
I Zw 18 stars compared to a local sample of evolved stars. 
Also stars that were not detected in F171M and F180M 
are included. Solid circles show IZwl8 intermediate-age 
stars, open circles show red (F110W — F160W > 0.15, 
and F110W — F205W > 1) young stars and open triangles 
show bluer young stars. The local sample is shown as solid 
triangles for oxygen-rich AGBs and supergiants, and open 
squares for carbon stars. 



in Sect. 3, only a small fraction of the stars are located 
near Ha sources, and we do not believe that this could be 
the main mechanism, but if it was, the net effect would be 
to make the intermediate-age population in IZwl8 older 
than it now appears. 

In figureOlwe compare the near-IR colours of stars in 
IZwl8, to a local sample of evolved stars from Langon 
and Wood IpOOO). A striking feature of the intermediate- 
age stars in I Zw 18 is that while they span a similar range 
in F160W-F205W, their F110W-F160W colours are bluer 
compared to the local AGB star sample, (see also Fig. 03 
and Sect. 14.21 for more details.). The observed difference 
between the stars in I Zw 18 and the local counterparts 
demonstrates that the low-metallicity stars with interme- 
diate mass and age, in IZwl8, have a F110W flux excess 
compared to local stars. On the other hand, young stars 
in I Zw 18 show bluer colours than the local sample for 
both colours, and their location in the F110W-F160W 
vs. F160W-F205W diagram is a continuation of the lo- 
cal oxygen-rich star sequence, reflecting the known effect 
of metallicity on stellar effective temperature and broad- 
band photometry. 
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4.2. Photometric properties of metal-poor stars of 
intermediate age and mass 

A caveat when modelling the near-IR CMDs of low- 
metallicity stellar systems, and for the identification of 
carbon star candidates, is the lack of representative spec- 
tra for metal-poor evolved stars in the stellar library 
used to transform the fundamental stellar parameters to 
observed quantities. The stellar library we have used is 
based on spectra of nearby stars, and is by construction 
biased toward metallicities consistent with the chemical 
evolution of our own Galaxy. The observed difference in 
F110W-F205W colours between stars of intermediate age 
and mass in I Zw 18 and their local counterparts may im- 
ply that the carbon star sequence in metal-poor systems 
would not be tightly constrained by a direct comparison 
with isochrones based on stellar libraries of local stars. 
In a metal-poor environment carbon stars are expected 
to appear at higher temperature compared to a metal- 
rich environment, and the progenitors of carbon stars will 
have higher effective temperature than their counterparts 
in the Large Magellanic Cloud (hereafter LMC) and the 
Milky Way. In addition, the convective dredge-up that en- 
riches the stellar atmosphere with carbon, is triggered ear- 
lier and is more efficient at low metallicity (see Mouhcine 
& Lancpn 2003 for a detailed discussion of these issues). 
However it is not an easy task to predict the distribution of 
carbon star colours. For Galactic and LMC carbon stars, 
JHK colours become redder as the C2 band strength, in- 
crease (Cohen et al. 1981). On the other hand, theoretical 
models show that the C2 band head is sensitive to the 
carbon-to-oxygen ratio (Gautshy-Loidl 2001). Synthetic 
evolutionary models for AGB stars predict that the fi- 
nal carbon-to-oxygen ratio increases as the initial stellar 
metallicity decreases (Mouhcine & Lancpn 2002, Marigo 
et al. 2002). Observationally, Matsuura et al. (2002) have 
shown that the equivalent widths of molecular absorption 
bands are systematically larger for carbon stars in the 
LMC than for carbon stars in the solar neighbourhood. 
They argue that the carbon-to-oxygen ratio is system- 
atically larger in LMC carbon stars. The extent of line 
blanketing, responsible for the distinction between car- 
bon stars and oxygen-rich stars, is related to the amount 
of carbon in the atmosphere. Thus, even if the progeni- 
tors of carbon stars are hotter at low metallicities, their 
spectrophotometric properties might be unaffected by the 
lower metallicity as long as the opacity in the atmosphere 
is dominated by molecules based on carbon and other el- 
ements whose abundances are not drastically affected by 
the low initial metallicity. 

Red F110W-F160W colours of local carbon stars are 
due to the effect of molecular blanketing, i.e., CN and 
C2 absorption bands. Atmospheric carbon abundance is 
larger at low metallicity, which may enhance the C2 and 
CN absorption bands. On the other hand, as a conse- 
quence of the extremely low abundance of nitrogen in 
IZwl8 ([N/O] = -0.7, [N/C] = -0.3, in solar units, 
Izotov and Thuan 1999), the abundance of molecules that 



contain nitrogen may be lower in the atmosphere of AGB 
stars in IZwl8 than in their local counterparts, lowering 
the CN-dominated opacity in F110W filter, but increasing 
the opacity in both the F160W and F205W filter wave- 
length intervals, where C2 dominates. This could give rise 
to the observed blue F110W-F160W and F110W-F205W 
colours. 

The occurrence of hot bottom burning in massive AGB 
stars, which convert a fraction of the newly dredged-up 
carbon into nitrogen via the CNO cycle, affects the abun- 
dance of both carbon and nitrogen in the atmosphere of 
such stars. It would be hard to derive any meaningful con- 
straints on the properties of the hot bottom burning at low 
metallicity on the basis of the small number of stars in our 
sample. Deeper observations of I Zw 18, and/or other low- 
metallicity systems, as well as detailed modelling of AGB 
stars, are needed to address this question. 

These findings illustrate that I Zw 18 constitutes an im- 
portant laboratory for understanding stellar physics and 
evolution at low metallicity. This is not only true for AGB 
stars - the discovery of Wolf-Rayet features in I Zw 18 a 
few years ago came as a surprise in view of the very low 
metallicity (Izotov et al. 119971 Legrand et al ll997l Brown 
et al. 12002) • 

4.3. Carbon star classification through medium-wide 
band photometry 

An alternative empirical approach to determine the nature 
of intermediate-age stars is to compare equivalent data for 
stars with known classification. In this subsection we try 
to classify the stars in I Zw 18 by comparing their photo- 
metric properties to those of a local sample from Lancpn & 
Wood (2000) that includes AGB stars (both oxygen-rich 
stars and carbon stars), red giants, and supergiants. 

The F180M-F171M index is a measure of the C 2 band 
strength for carbon stars or the H2O band strength for 
oxygen rich stars. It is able to separate between evolved 
AGB stars and other types of evolved stars, but it cannot 
separate the two sub-populations of AGB stars (Alvarez 
et al. 2000). Fortunately, combining this index with other 
information can do a good job in separating both AGB 
sub-populations. Indeed, local carbon stars split from the 
rest of the local sample in the F180M-F171M vs. F171M- 
F160W diagram. The F171M-F160W index measures the 
slope of the spectrum between the centre and the edge of 
the F160W filter. It has a numerically larger for oxygen- 
rich stars than for carbon stars as a consequence of the 
presence of broad H2O absorption band, while for car- 
bon stars the flux at the edge of the filter is larger than 
in the centre. As the spectral features of carbon stars 
are not strongly sensitive to stellar temperature (Lancpn 
& Mouhcine 2002), the behaviour of carbon stars in the 
(F180M-F171M) vs. (F171M-F160W) diagram should in 
principle be independent of metallicity. To identify car- 
bon star candidates, we consider stars that are classified 
as intermediate-age stars in the near-IR CMDs (Fig. 2), 
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Fig. 4. Diagnostic diagrams for IZwl8 stars with available F171M and F180M photometry. Upper left panel shows 
F180M-F171M vs. F171M-F160W compared to a local sample of evolved stars: triangles show spectroscopically 
classified carbon stars, and circles show oxygen rich stars, giant stars, and supergiants. Filled squares show intermediate- 
age stars. The arrow indicates that one of these with F171M — F160W = —1.5 is located to the left of the border of 
the plot. The upper right panel shows F180M-F171M vs. F110W-F205W diagram. The bottom left panel displays 
F205W vs. F180M-F171M diagram. The bottom right panel shows F110W-F160W vs. F160W-F205W diagram. Stars 
classified as young stars are shown as open squares. For the last three panels, carbon star candidates are shown as 
filled squares, an oxygen-rich intermediate-age star is shown as a filled circle. 

i.e., older than 100 Myr. Six of these have available F180M the Langon & Wood (2000) local sample, where triangles 
and F171M photometry. show carbon stars, and circles show all other oxygen rich 

In the upper left panel of Fig.H filled squares show stellar types, i.e. giants, supergiants, and oxygen-rich AGB 
the distribution of intermediate-age IZwl8 stars in the stars. As the diagram shows, the combination of F180M- 
F180M-F171M vs. F171M-F160W diagram, compared to F171M and F171M-F160W indices can efficiently separate 
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Fig. 5. Similar to Fig.[21 but restricted to stars with available medium-wide band photometry. Filled squares show 
carbon star candidates. The filled circle shows an oxygen-rich intermediate-age star (see text for more details). 



spectroscopically identified oxygen-rich and carbon stars, 
as the carbon stars have systematically smaller F171M- 
F160W at a given F180M-F171M. The six intermediate- 
age stars in I Zw 18 fall into two groups in this diagram: 
Five stars, generally with positive F180M-F171M index, 
fall close to, or to the left of, the local carbon stars. Only 
one intermediate-age star in I Zw 18 falls close to the lo- 
cation occupied by local supergiant and oxygen-rich AGB 
stars. The other I Zw 18 stars in this region have too blue 
broad band colours and/or are too bright to be classified 
as intermediate-age stars, although there are two border- 
line cases. Bright stars with F110W - F205W < 1 tend 
to have F171W — F160W ~ 0, while fainter stars span 
a larger range in F171W-F160W. There seems to be a 
anti-correlation between the F180M-F171M and F171M- 
F160W indices for intermediate-age IZwl8 stars: the for- 
mer index gets stronger on average as the latter index get 
weaker, and the same trend is seen for local/LMC carbon 
stars. Hence, according to the upper left panel of Fig.^J 
the stars in I Zw 18 that are located in the region occu- 
pied by local carbon stars are likely carbon stars. The 
single intermediate-age star located close to oxygen-rich 
stars is likely an oxygen-rich AGB star. Below, we use 
additional diagnostic diagrams to confirm this classifica- 
tion. Note that, as discussed in Sect. 3, this classification 
is completely insensitive to any possible contamination of 
nebular emission. 

The upper right panel of Figure^] shows the 
F180M-F171M vs. F110W-F205W diagram for stars in 
IZwl8. The figure shows again the dichotomy between 
intermediate-age/mass stars and young/massive stars. For 
stars redder than F110W — F205W « 1, intermediate- 



mass stars tend to have positive F180M-F171M index, 
while supergiant stars have in average negative F180M- 
F171M index. Only two supergiants (with Mp205W < — 9) 
have a positive F180M-F171M index, including the most 
luminous star in the sample, and of all the young stars, 
only 6 out of 22 have F180M - F171M > 0. The fig- 
ure shows also the inability of the F180M-F171M index 
alone, to separate between oxygen-rich and carbon stars. 
The 180M-171M index, as expected for intermediate-mass 
stars, gets stronger on average as the stars get fainter as 
shown in the bottom left panel of Fig.Q] 

In the bottom right panel of Fig.QJ we show that car- 
bon star candidates tend to have redder F160W-F205W 
colour at a given F110W-F160W colour, similar to what is 
seen for LMC/Galactic carbon stars (Fig. - In principle, 
nebular contamination would give a similar behaviour, but 
this explanation is unlikely since the stars with the most 
deviant colours are the carbon star candidates which are 
not located near any Ha sources. This rather supports 
the molecular blanketing interpretation, as opposed to a 
simple effective temperature effect. Other stars that show 
slightly redder F160W-F205W colours at a given F110W- 
F160W, but not classified as carbon stars, are too blue 
and/or too luminous to be intermediate-mass stars. Fig.03 
shows the same CMDs as in Fig.|2 but restricted to stars 
that have been detected in medium- wide filters. It is clear 
that almost all stars older than ~ 100 Myr are carbon 
stars. This is consistent with the theoretical models that 
predict that carbon stars dominate intermediate-age stel- 
lar populations at low metallicities. 

We have shown that the F171M-F160W vs. F180M- 
F171M diagram can separate probable carbon stars from 
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Fig. 6. Spatial location of stars in I Zw 18. Left: co-added F110W, F160W and F205W frame with logarithmic intensity 
scaling. The orientation is indicated by the compass in the lower right corner, and the length of the arrow is 2". Right: 
Filled (red) circles show intermediate-age stars detected in F110W, F160W and F205W, where those brighter than 
F205W = 22 have slightly larger symbols. Those intermediate-age stars with available F171M and F180M photometry 
that have been identified as carbon star candidates have been highlighted by rays emanating from them. The candidate 
oxygen-rich star is marked with a plus (and shown in blue). Open circles show young (< 100 Myr) stars, where stars 
brighter than F205W= 22 are shown by larger symbols. Stars detected in F110W and F160W only are shown as small 
crosses, if F110W-F160W< 0.7, or as small dots otherwise. One source suspected to be a background galaxy is in 
addition labelled with a "G" . 



other late type stars. If we relax the requirement of hav- 
ing photometry in both the F171M and F180M filters, 
and look also at stars accurately detected in F171M only, 
the sample is increased with four more stars. Two of 
these are faint and blue and most probably young stars, 
whereas two have broad band photometry consistent with 
intermediate-age stars. The two latter stars have F171M- 
F160W = —0.1 and —0.3, hence most consistent with local 
O-rich stars. However, without F180M-F171M photome- 
try, this classification becomes tentative. Whereas our five 
carbon star candidates separate well from the rest of the 
stars in a F171M-F160W vs. F160W-F205W diagram this 
is not true for carbon stars in general as can be understood 
from the upper left of Fig. 0] 

4.4. Spatial location of stars in IZwl8 

The spatial location of both young supergiant stars and 
intermediate-age stars are shown in Fig. Open circles 
show young stars, and filled circles show intermediate- age 
stars. Carbon stars candidates and the candidate AGB 



oxygen-rich star are highlighted. The figure shows that 
luminous supergiants are concentrated to the NW and SE 
regions, where the star formation is active, while AGB 
stars are more uniformly distributed. 

4.5. Pixel photometry of unresolved stars 

As described in Sect.[3]we have analysed the colours of in- 
dividual pixels in the SE region. The result for a 30 arcsec 2 
area is shown in Fig. (noisy pixels and pixels contam- 
inated by detected stars have been excluded, for details 
see Sect. [3J, together with the photometry of the stars 
discussed above. The region is approximately quadratic 
with a centre at (x, y) = (—2, 0) in the coordinate system 
of Fig. EJ 

We have seen that the F171M-F160W vs F180M- 
F171M colour-colour diagram is an efficient tool for dis- 
criminating carbon stars and other class of objects. Pixels 
that occupy the location of spectroscopically identified 
carbon stars in this diagram are then considered to be 
dominated by carbon stars. Figure 7 shows that individ- 
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Fig. 7. Pixel photometry of south-east region Left: F171M-F160W vs F180M-F171M photometry. Each dot shows 
the photometry of individual pixels (binned to the original 0.075 arcsec/pixel scale of NIC2) in the south-east region 
with a luminosity in F171M and F180M brighter than 26.6 magnitudes, and excluding the point sources detected these 
filters. The filled squares show the corresponding photometry for stars in IZwl8 (see the upper left panel of Fig.0J). 
Right: F160W-F205W vs F110W-F160W photometry, for the same pixels as presented in the left plot. The filled 
squares show all stars in I Zw 18 (see Fig. [3} 



ual pixels have a broad range of colours indicating that 
the unresolved stellar population contains both young and 
intermediate-age stars. Interestingly, the carbon star re- 
gion (see upper left of Fig. 4) is well populated. The in- 
tegrated colours for pixels in the carbon star region are 
F171M-F160W = -1.1 and F180M-F171M= 0.0, and 
have a total F160W magnitude of 21.6. For single stars de- 
tected in all 5 passbands, the limiting magnitude is F160W 
?» 23.5 meaning that the unresolved pixels could corre- 
spond to, on the order of 10 carbon stars with F160W~ 24, 
or Mpi6ow = -6.5. 

Taking the area of I Zw 18 to be 14" x 8" (the angular 
size at = 25 mag/arcsec 2 , Papaderos et al. I2002|l . the 
surface density of detected carbon stars in IZwl8 as a 
whole is ~ 15 per kpc 2 . The SE region has an area which 
is 27% of the total assumed area and contains one of the 5 
identified carbon stars. If the SE region is assumed to be 
representative for the total AGB population of I Zw 18, the 
pixel photometry of this region suggests that the surface 
density may be on the order of 100 per kpc 2 , and the total 
number of carbon stars on the order of 50, but of course 
this sensitively depends on the assumed luminosity func- 
tion of carbon stars. This would mean that IZwl8 has a 
carbon star surface density rather typical for faint metal- 
poor dwarf galaxies (see Groenewegen 1999). The absolute 
magnitude of I Zw 18 is My as —14, but this is boosted by 
the presence of many young stars. Assuming that the av- 



erage metallicity of I Zw 18 stellar content is similar to 
its gas phase metallicity, i.e., Z = 0.02 Zq, IZwl8 falls 
on the correlations of the total number of carbon stars 
against galaxy luminosity, and the normalised number of 
carbon stars to galaxy luminosity against galaxy metal- 
licity respectively, as defined by Local Group galaxies 
(Groenewegen ^23- Consequently carbon star content of 
I Zw 18 appears to be similar to these observed for metal- 
poor Local Group dwarf galaxies. Given the significant dif- 
ference between the luminosity- weighted ages of the Local 
Group dwarf galaxies and IZwl8, this may indicate that 
the scaling correlations that regulate the properties of car- 
bon star content in galaxies could be metallicity rather 
than age. 



5. Conclusions 

We have presented new broad and medium-wide band 
HST/NICMOS imaging of the very metal-poor blue 
compact galaxy IZwl8. The data indicate that an 
intermediate-age stellar component is indeed present, and 
that the galaxy was not formed during the last 100 Myr. 
We emphasise the usefulness of combined broad- and 
medium-wide band near-IR photometry for detecting and 
identifying evolved intermediate-age stars whose energy 
distribution peak in the near-IR. 
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Based on comparison with isochrones and on empir- 
ical matching of stars in IZwl8 to morphological fea- 
tures of known stellar populations in different diagnos- 
tic diagrams we have identified carbon star candidates. 
Out of the stars classified as intermediate-age through 
F110W (~J), F160W (~H), and F205W (~K) photom- 
etry, six have available medium- wide band photometry in 
the F171M and F180M filters, and five of these are prob- 
able carbon stars. This suggests that the bright AGB star 
population in I Zw 18 is dominated by carbon stars. At a 
distance of more than 10 Mpc, these are the most distant 
resolved carbon stars discovered as yet. Investigating the 
broad and medium- wide band photometry for pixels in the 
south-east region we find support for more carbon stars 
lurking behind the single star detection threshold. 

We find that intermediate-age stars in IZwl8 show 
relatively blue F110W-F160W colours compared to local 
AGB stars. The opacity source in AGB stars for the wave- 
length region covered by the F110W filter is dominated by 
the CN molecule, and we speculate that the blue colours 
may be due to the very low nitrogen abundance in I Zw 18. 
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